Nonlinear Supression of Tunneling and Strong Localization in Bose-Einstein 
Condensates with Spatially Inhomogeneous Interactions 
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We study the properties of the ground state of Bose-Einstein condensates with spatially inhomo- 
geneous interactions and show that the atom density experiences a strong localization at the spatial 
region where the scattering length is close to zero while tunneling to regions with positive values of 
the scattering length is strongly supressed by the nonlinear interactions. 
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The experimental generation of Bose-Einstein conden- 
sates (BEC) with ultracold dilute atomic vapors [l| has 
turned out to be of exceptional importance for physics. 
The formation of a condensate occurs when the temper- 
ature is low enough and most of the atoms occupy the 
ground state of the system. This process is visible both in 
momentum and in real space due to the spatial inhomo- 
geneities exhibited by the order parameter on a macro- 
scopic scale because of the trapping potentials. 

The properties of the ground state of trapped BECs 
are well known. In the mean field limit simple analyt- 
ical expressions are available in the Thomas-Fermi ap- 
proximation [2] and beyond Q. These approximations 
and direct numerical simulations describe accurately the 
properties of the experimentally found ground states Q . 

Nonlinear interactions between atoms in a Bose- 
Einstein condensate are dominated by the two-body col- 
lisions that can be controlled by the so-called Feschbach 
resonance (FR) management [5j. The control in time of 
the scattering length has been used to generate bright 
solitons @ and induce collapse Q and has been the ba- 
sis for theoretical proposals to obtain different nonlinear 
waves [H, 0, [Toj| . Moreover, interactions can be made spa- 
tially dependent by acting on either the magnetic field or 
the laser intensity (in the case of optical control of FRs 
[Tlj ) acting on the Feschbach resonances. This possibil- 
ity has motivated many theoretical studies on the behav- 
ior of solitons in Bose-Einstein condensates (BECs) with 
spatially inhomogeneous interactions [l2|, [jjj [li], H, G3 • 

In this paper we study the effect of specific spatially in- 
homogeneous interactions on the ground state of a BEC 
in the mean field limit. We will show that when the 
scattering length is non- negative, a striking localization 
phenomenon of the atom density occurs at the regions 
where the scattering length vanishes. By tuning appro- 
priately the control (magnetic or optical) fields this phe- 
nomenon can be used to design regions with large particle 
densities and prescribed geometries. Another interesting 
phenomenon to be studied in this letter is the nonlinear 
limitation of tunneling of atoms to the regions in which 
the interactions are stronger. 



The model and its Thomas-Fermi limit.- The ground 
state of a BEC in the mean field limit is the real, positive 
solution of the Gross-Pitaevskii equation 



\ ( j ) = --A^ + V{x)4> + g{x)\4>\ 2 4>, (1) 



which minimizes the energy 



E{<f>) = [ 

JR 3 



(2) 



under the constraint of a fixed number of particles N — 
J R3 |(/>| 2 . Eq. (pp) is written in nondimensional units where 
the coordinates x and time t are measured in units of 
ao = y^h/muj and 1/cj, respectively, while the energies 
and frequencies are measured in units of Huj and uj re- 
spectively, ou being a characteristic frequency of the po- 
tential. Finally, g(x) = 47ra(x)/ao is proportional to the 
local value of the s-wave scattering length a. 

The Thomas-Fermi approximation proceeds by ne- 
glecting the kinetic energy or equivalently the term pro- 
portional to A(j) in Eq. (pQ). In the case of spatially ho- 
mogeneous interactions g(x) = go this leads to </>tf(#) = 
[(A — V(x)) /go] 1 ^ 2 - When the interactions are spatially 
inhomogeneous the same formal manipulation leads to 



4> TF {x) = y/(X-V(x))/g(x), 



(3) 



which diverges on the set G = {x G M 3 : g(x) = 0}. 
Obviously, in the vicinity of G the Thomas-Fermi ap- 
proximation is not correct but anyway this divergence is 
a first indication of a phenomenon to be studied in this 
paper: the tendency of the ground state in Bose-Einstein 
condensates with spatially inhomogeneous interactions to 
localize strongly on the regions where the scattering length 
is close to zero provided the system is sufficiently nonlin- 
ear, i.e. for sufficiently large values of gN . 

A simple example.- Let us first consider an exactly 
solvable "toy" example which displays the main features 
of the phenomena to be studied in this paper: a quasi- 
one dimensional BEC in a box, i.e. setting V(x) = 0, 
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FIG. 1: [Color online] Spatial distribution of the density 
n(x) — \(j)(x)\ 2 for the ground state solutions of Eq. Q 
with L = 2, a = 1 and different values of the chemical po- 
tential A corresponding to different values of N (a) From the 
lower to the upper curve goN = 2,3.55 (corresponding to 
A*) and goN = 6. (b) From the lower to the upper curve 
g N = 25, 50, 100, 200. 



4>(x = ±L) = with scattering length given by g(x) = go, 
for \x\ < a and g(x) = 0, for \x\ > a. 
In this simple case, Eq. (pp) becomes 



<j> xx + 2\(j> = 2# </> 3 , \x\<a, 
<j> xx + 2\(j) = 0, a < \x\ < L, 



(4a) 
(4b) 



and its positive, even solution, satisfying the boundary 
conditions 0(=LL) = 0, </>'(0) = is given by 



<t>(x) =Csin V2X(x-L) 



, a < \x\ < L, 



(5) 



while for \x\ < a the solutions are given by 

^ sn (xV\a + 5; fc 2 ) , A < A* 
4>{x) = { 7r/(2 v / 2^|o-L|) , A = A*, 

^ dc (xV\^; fc 2 ) , A > A* 



(6) 



where a(k) = 2/(1 + & 2 ), sn and dc are two of the stan- 
dard Jacobi elliptic functions and k is the elliptic mod- 
ulus. Both the elliptic modulus and amplitude C can 
be obtained from the matching conditions for </>(a) and 
4>'{a). These conditions also give A* = 7r 2 / [8 (a — L) 2 ]. 
Finally, the cutoff value of the chemical potential A c can 
be obtained from the condition of maximum slope at 
x = a which leads to A c = 7r 2 /[2(a — L) 2 ]. In that case 
the number of particles in the condensate is infinite, since 
the amplitude in the outer region C — > oo. 

The spatial profiles of the ground state density for dif- 
ferent values of go^V" shown in Fig. [T] support our con- 
jecture based on the Thomas-Fermi solution, i.e. the 
existence of a strong localization of the atom density in 
the region where the interactions vanish. 

It is also remarkable that the atom density in the in- 
ner part of the domain, i.e. the region where there are 
nonlinear interactions, remains almost constant indepen- 
dently on the number of particles in the condensate once 
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FIG. 2: [Color online] (a) Ratio between the maximum atom 
density and the atom density at x = as a function of the 
scaled scattering length g* on the spatial region |x| < a (for 
go = 1) (b-c) Scaled atom density profiles for (b) g* = and 
(c) g* =0.2 in both cases for a total scaled number of particles 
N = 1000. The insets show the profile of g(x). 



a certain critical density is achieved. This region corre- 
sponds to the nonlinear analogue of the classically for- 
bidden region in ordinary potentials and is energetically 
less favourable due to the extra repulsive energy provided 
by the nonlinear interactions. However, the tunneling of 
atoms in this region is essentially limited to a constant 
value, independently of the number of atoms, which dif- 
fers essentially from ordinary tunneling. 

The supression of tunneling depends strongly on the 
value of the scattering length in the inner region, that 
we have taken to be zero up to now. If instead we set 
g(x) = g* when \x\ < a and study the dependence of the 
ratio between the maximum atom density and the atom 
density at x = (which is a measure of the amplitude 
of the tunneling), we find a strong dependence on this 
parameter as shown in Fig. EJa). 

This effect is also seen in the atom density profiles 
when comparing the cases with g* = [Fig. [5(b)] and 
g* = 0.2 [Fig. [2(c)] for N = 1000. Larger values of TV 
lead to a stronger effect. 

Numerical results.- Let us now consider Eq. (pQ) in 
more realistic quasi one-dimensional scenarios by includ- 
ing a potential V(x) = 0.02x 2 . We have computed 
numerically the ground state for scattering lengths of 
the form go(x) = l : gi(x) = exp(— x 2 /200), g2(x) = 
exp(— x 2 /50). These choices allow us to study different 
degrees of localization of the interactions starting from 
the case of no localization. Our results are summarized 
in Fig. El 

In Fig. [3(a) we observe how the maximum density 
(n(x) = \(j)(x)\ 2 ) increases drastically for spatially decay- 
ing nonlinearities (blue and red curves) as a function of 
the number of effective number of particles in the quasi- 
one dimensional condensate N. This amplitude growth is 
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FIG. 3: [Color online] Ground states of Eq. JTJ for V(x) = 
0.02x 2 and go(x) = l,gi(x) = exp(— x 2 /200), and g2(x) — 
exp(— x 2 /50), for different values of the scaled number of par- 
ticles N . (a) Maximum particle density max x n(x) and (b) 
width W 2 — f x 2 n(x)dx/N for go (green), gi(x) (red) and 
92(x) (blue), (c) Spatial profiles of y/n(x) for N = 10 (blue), 
N = 100 (green), N = 1000 (red) for g(x) = gi(x). The 
dashed black line is the ground state with homogeneous in- 
teractions and N = 1000. (d) Same as (c) but for N = 10000 
(blue), and N = 40000 (green), in comparison with the case 
of spatially homogeneous interactions (dashed black lines). 
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FIG. 4: [Color online] Ground state maximum particle density 
(a) and width W = fx 2 n(x)dx/N (b) for V(x) = 0.02x 2 and 
go(x) — 1 (red lines) and g±(x) — (1 — 0.001x 2 )+ (blue lines). 
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FIG. 5: [Color online] Ground states for V(x) — 0.02x 2 and 
g(x) = 1 — exp(— x 2 /50) (shown in panel (c)). (a) y/n(x) for 
N = 10 (blue line), N = 100 (green line) and N = 1000 (red 
line), (b) Condensate width as a function of N. 



due to a strong localization effect near the region where 
g(x) vanishes as it is seen in Fig. Etc,d). In contrast, 
the condensate density for spatially homogeneous inter- 
actions grows slowly according to the Thomas- Fermi pre- 
diction maxn(x) ex TV 2 / 3 . When the number of particles 
is small, the size of the atomic cloud is smaller than the 
localization region of g(x). For larger values of TV the 
ground state extends beyond the localization region of 
g(x) and the atom density becomes more and more lo- 
calized near its edge. This effect is more clear for larger 
number of particles and is accompanied by a saturation 
in the amplitude growth in the region where g(x) is far 
from zero [Fig. 13(d)]. 

In the case of spatially homogeneous interactions the 
width grows according to the law W oc TV 1 / 3 [Fig. [3(b)]. 
When interactions are spatially dependent and due to the 
localization of the amplitude close to the zero of g(x) the 
width growth saturates for large values of TV to a value 
depending on the size of g(x). 

These effects are even more clear when the nonlinear- 
ity decays to zero faster. For instance, taking a nonlinear 
coefficient given by g±(x) = (1 — 0.001x 2 )+ (i.e. an in- 
verted parabola with maximum amplitude g = 1 at x = 
and zero values for \x\ > 31.6) as shown in Fig. [H we 
see that the maximum densities are even higher than for 
<7i 0*0 5 #2 0*0 • In this case, in comparison with our first 
simple example, we can see an even stronger localization 
since the existence of the potential makes energetically 



more favourable the localization close to the point where 
the interactions vanish. 

Finally, in Fig. [5] we get again the localization phe- 
nomenon but now near x = for localized interactions 
given by g(x) = 1— exp(— x 2 /50) and V(x) = 0.02x 2 . It is 
interesting how localized the density becomes to "avoid" 
penetrating into the regions with appreciable values of 
the scattering length. 

Rigorous results.- In Ref. [17] the equation 

— Au = Xu — a(x) (u r + /(x, u)) u, (7) 

was considered on a bounded domain Q C M N with 
C 1 ' 1 boundary under the conditions u\on = 0. r > 
0, a(x) G C(Q) is nonnegative and the open set 
D = {xgO: a(x) > 0} satisfies D C ft and pos- 
sesses a finite number of C 1 connected components 
Dj,l < j < I such that A H Dj = if i / j. 
The function / : Cl x [0, +oo) — > R, satisfies /, f u = 
df/du G C (Q, x [0, oo);R) and the growth conditions 
littMo /('j u ) = °> lim nToo f u (-, u)/^- 1 = 0, uniformly in 
Q. In the case when / = 0, Eq. ([7j) is the logistic model of 
population dynamics [18[ where Q is the region inhabited 
by the species u, A measures its birth rate and a(x) mea- 
sures the capacity of Q to support the species u. Given an 
open subset Q± C O with a finite number of components 
let us denote by crp 1 [—A] the minimum of the principal 
eigenvalues of the operator —A on each of the compo- 
nents of Qi subject to homogeneous Dirichlet boundary 
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conditions on dQ i it was proven [ItJ that Eq. ([7j) has an 
unique positive solution for A G (crp[— A], crp° [— A]) and 
that no such solution exists for A > <rp°[— A]. Moreover, 
it grows to infinity when A | [—A] on compact subsets 
of and stabilize in D to the minimal weak solution of 
Eq. (O with singular boundary conditions u\qd = oo. 

From Eq. (JJ) we can obtain our quantum mechani- 
cal problem described by Eq. ([Tj) by choosing f(x,u) = 
0, a(x) = 2g(x),r = 2. However, the results of Ref. [17[ 
essentially require a bounded domain and a continuous 
g(x) and thus are not directly applicable to the examples 
discussed in this paper. However, they support the local- 
ization phenomenon and the saturation of the tunneling 
amplitude in formally related problems. 

Conclusions and discussion.- We have studied the 
ground state of Bose-Einstein condensates with repulsive 
spatially varying interactions and shown that the atom 
density localizes dramatically near the regions where the 
interactions are close to zero when the nonlinear interac- 
tions are relevant enough. This behavior does not depend 
on the spatial dimensionality and can be used to design 
very effectively one, two or three-dimensional spatial dis- 
tributions with large values of the atom density by acting 
on the control field, e.g. using micromagnets to change 
locally the magnetic field or localized laser beams. This 
behavior is not so easy to achieve with ordinary poten- 
tials acting on the condensate. Our results could also be 
useful for atom lithography, atom beam guiding or other 
applications where a precise control of the positioning of 
large values of the density of an atomic cloud is required. 

We have also found a limitation of tunneling to the 
regions with larger values of the scattering length which 
is independent on the number of particles and differs es- 
sentially from classical tunneling in ordinary potentials. 

These phenomena may appear in multicomponent con- 
densates that offer a much wider range of possibilities for 
their ground states depending on the type of components 
and their relative interactions [l9(] . It would be interest- 
ing to study the situation in which the variations of the 
scattering length make the components to be miscible 
or inmiscible in different spatial regions and what is the 
geometry of the resulting domains. 

These phenomena could also be observed in nonlin- 
ear optical systems. Although ordinary optical materials 
have small nonlinear it ies, there are media with a giant 
nonlinear response [20[ as it happens when a probe laser 
beam propagates in a medium with transparency induced 
electromagnetically by a second coupling field. There- 
fore, managing the parameters it could be posible to find 
an optical version of the phenomena presented here. 
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